#/usr/bin/Rscript

on_100_20 <- read.table("on/profits.dat")
op_100_20 <- read.table("rf/profits.dat")
ut_100_20 <- read.table("ut/profits.dat")
print("100 20 on")
shapiro.test(on_100_20$V1)
print("100 op")
shapiro.test(op_100_20$V1)
print("100 ut")
shapiro.test(ut_100_20$V1)

#QQplots

jpeg("normalidade_100.jpg")
par(mfrow=c(1, 3))

qqnorm(on_100_20$V1, main="on 100 10")
qqline(on_100_20$V1)
qqnorm(op_100_20$V1, main="op 100 10")
qqline(op_100_20$V1)
qqnorm(ut_100_20$V1, main="ut 100 10")
qqline(ut_100_20$V1)

jpeg("homocedasticidade_100.jpg")
par(mfrow=c(1, 3))

plot(on_100_20$V1, main="on 100 10")
plot(op_100_20$V1, main="op 100 10")
plot(ut_100_20$V1, main="ut 100 10")

dev.off()

